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We consider the production of W and Z bosons in hadron collisions. We present 
a selection of numerical results obtained through a fully exclusive calculation up to 
next-to-next-to-leading order (NNLO) in QCD perturbation theory. We include the 
7-Z interference, finite-width effects, the leptonic decay of the vector bosons and the 
corresponding spin correlations. The calculation is completely realistic, since it allows 
us to apply arbitrary kinematical cuts on the final-state leptons and the associated 
partons, and to compute the corresponding distributions in the form of bin histograms. 

The production of Icpton pairs through the Drell-Yan (DY) mechanism [1] was the first 
process where parton model ideas developed for deep inelastic scattering were applied to 
hadronic collisions, and lead to the discovery of W and Z bosons [HIS]. 

It is thus not surprizing that the production of vector bosons is central in physics studies 
at hadron colliders. These processes have large production rates and clean experimental 
signatures, given the presence of at least one high-pT lepton in the final state. Studies of 
the production of W bosons at the Tevatron lead to precise determinations of the W mass 
and width [5]. Vector boson production is also expected to provide standard candles for 
detector calibration during the first stage of the LHC running. 

For these reasons it is important to have accurate theoretical predictions for the vector- 
boson production cross sections and the associated distributions, and such a task requires 
detailed computations of radiative corrections. The QCD corrections to the total cross 
section [5] and to the rapidity distribution [6l [7] of the vector boson are known up to the 
next-to-next-to-leading order (NNLO) in the strong coupling as. The fully exclusive NNLO 
calculation, including the leptonic decay of the vector boson, has been completed more 
recently [HI [5]. Full electroweak corrections at 0{a) have been computed for both W [TU] 
and Z production [TT] . 

In this contribution we discuss a recent computation of the NNLO QCD corrections to 
vector boson production in hadron collisions [S]. 

The evaluation of higher-order QCD corrections to hard-scattering processes is well 
known to be a hard task. The presence of infrared (IR) singularities at intermediate stages 
of the calculation prevents a straightforward implementation of numerical techniques. In 
particular, NNLO differential calculations are a rarity due to their substantial technical 
complications. In e+e" collisions, NNLO differential cross sections are known only for 
two [I2I [13] and three jet production [HI [15]. At hadron colliders fully differential cross- 
sections have been computed only for Higgs production in gluon fusion [TU [171 El [H] j and 
the Drell-Yan process [SI [9] . It is interesting to note that the amplitudes relevant for vector 
boson production at NNLO have been known for at least 15 years |S] before the first fully 
exclusive computation could be completed [8]. 
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Figure 1: Rapidity distribution of the bo- 
son at the Tevatron. The NNLO result (blue) 
is compared to the NLO band (red) and to the 
NNLO prediction of Ref. [7]. 



The calculation [5] we discuss here is 
based on an extension of the subtraction 
formahsm [20l[2T] to NNLO that can be ap- 
plied to the production of colourless high- 
mass system in hadron collisions |17j . The 
calculation parallels the one recently com- 
pleted for Higgs boson production [T71 [TH] , 
and it is performed by using the same 
method. We include the j-Z interference, 
finite- width effects, the leptonic decay of the 
vector bosons and the corresponding spin 
correlations. 

In the following we present some numer- 
ical results for W and Z production at Teva- 
tron energies. We consider np = 5 massless 
quarks in the initial state, and, in the case 
of production, we use the (unitarity 
constrained) CKM matrix elements Vud = 
0.97419, Vus = 0.2257, Kb = 0.00359, 
Vcd = 0.2256, Vcs = 0.97334, Kb = 0.0415 from the PDG 



2008 [22]. In the case of Z 
production, additional Fcynman diagrams with fermionic triangles should be taken into ac- 
count. Their contribution cancels out for each isospin multiplet when massless quarks are 
considered. The effect of a finite top-quark mass in the third generation has been considered 
and found extremely small [26], so it is neglected here. 

As for the electroweak couplings, we use the so called scheme, where the input 
parameters are Gp , rnz, niw- In particular we use the values Gp = 1.16637 x 10~^ GeV~^, 
mz = 91.1876 GcV, Tz = 2.4952 GeV, mw = 80.398 GcV and Tw = 2.141 GcV. We 
use the MSTW2008 [23] sets of parton distributions, with densities and as evaluated at 
each corresponding order (i.e., we use {n + l)-loop as at N"LO, with n = 0,1,2). The 
renormalization and factorization scales are fixed to the value = fip = my, where my 
is the mass of the vector boson. 

We start the presentation of our results 
by considering the production of an on-shell 
boson at the Tevatron. When no cuts 
are applied our numerical program allows 
an independent computation of the rapidity 
distribution of a vector boson up to NNLO 
[7]. To compare with Ref. [7], in Fig. [1] we 
show the rapidity distribution of the 
obtained by using the MRST2001 partons 
[11I2S]. The blue histogram is the NNLO 
prediction; in red we also show the NLO 
band, obtained by varying ^p ~ be- 
tween l/2mw and 2mw. The solid curve 
is the (scaled) NNLO prediction extracted 
from Fig. 10 of Ref. [7]. The two NNLO 
results appear to be in good agreement. 
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Figure 2: Distributions in p^n 
the Z signal at the Tevatron. 
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We next consider the production of e+e^ pairs from Z/j* bosons at the Tcvatron. For 
each event, we classify the lepton transverse momenta according to their minimum and 
maximum values, pTmin and pTmax- The leptons are required to have a minimum px of 
20 GeV and pseudorapidity \ri\ < 2. Their invariant mass is required to be in the range 
70 GeV < TTie+e- < 110 GeV. The accepted cross sections are (Jlo = 103.37 ± 0.04 pb, 
aNLO = 140.43 ± 0.07 pb and aNNLO = 143.86 ± 0.12 pb. 

In Fig. [2] we plot the distributions in 
PTmin and pTmax at LO, NLO and NNLO. 
We note that at LO the pTmin and pTmax 
distributions are kinematically bounded by 
PT < <3max/2, where Qmax = 110 GeV is 
the maximum allowed invariant mass of the 
e~^e~ pairs. The NNLO corrections have a 
visible impact on the shape of the pTmin and 
PTmax distribution and make the prmin dis- 
tribution softer, and the prmax distribution 
harder. 

We finally consider the production of a 
charged lepton plus missing px through the 
decay of a boson {W = , W~) at the 
Tevatron. The charged lepton is required to 
have pt > 20 GeV and |?7| < 2 and the miss- 
ing pt of the event should be larger than 

25 GeV. We define the transverse mass of the event as tot = \J 2plpp^^^^{l — cos 4>) , where (f> 
is the angle between the the pr of the lepton and the missing px- The accepted cross sections 
are aLO = 1-161 ± 0.001 nb, aNLO = 1-550 ± 0.001 nb and aNNLO = 1-586 ± 0.002 nb. In 
Fig. [3] we show the tot distribution at LO, NLO and NNLO. We note that at LO the distri- 
bution has a kinematical boundary at tot = 50 GeV. This is due to the fact that at LO the 
W is produced with zero transverse momentum: therefore, the requirement p^^^^ > 25 GeV 
sets Tot > 50 GeV. Around the region where tot = 50 GcV there are perturbativc instabil- 
ities in going from LO to NLO and to NNLO. The origin of these perturbativc instabilities 
is well known |27| : since the LO spectrum is kinematically bounded by tot > 50 GeV, each 
higher-order perturbativc contribution produces (integrable) logarithmic singularities in the 
vicinity of the boundary. We also note that, below the boundary, the NNLO corrections to 
the NLO result are large; for example, they are about -1-40% at tot ~ 30 GeV. This is not 
unexpected, since in this region of transverse masses, the 0{as) result corresponds to the 
calculation at the first perturbativc order and, therefore, our 0{ag) result is actually only 
a calculation at the NLO level of perturbativc accuracy. 

We have discussed some selected results of a calculation of W and Z boson production 
up to NNLO in QCD perturbation theory. The calculation is implemented in a parton level 
event generator and it is particularly suitable for practical applications to the computation 
of distributions in the form of bin histograms. A public version of our program will be made 
available in the near future. 



Figure 3: Distributions in pTn 
the Z signal at the Tevatron. 
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